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I.  Introduction 


High  quality  topographic  models  and  maps  of  land  cover  are  critical  for  a  wide  range  of 
applications  including  strategic  military  operations,  disaster  response,  and  environmental 
monitoring.  These  requirement  are  evidenced  by  the  NASA  SRTM  ICESat,  and  EO-1 
missions,  as  well  as  advances  in  airborne  synthetic  aperture  radar  (SAR)  interferometry,  light 
detection  and  ranging  (LIDAR),  and  improved  multispectral  and  hyperspectral  sensors.  The 
goal  of  this  research  project  was  to  develop,  implement,  and  evaluate  computational  methods 
for  extracting  information  from  multiple  sources  of  remotely  sensed  data  for  thematic  and 
topographic  mapping  of  geologic  features.  The  research  focused  on  three  primary  areas: 

•  Integration  of  multisource  topographic  information  in  a  multiresolution 
framework 

•  Development  of  supervised  classification  methods  for  multispectral  and 
hyperspectral  data. 

•  Extraction  of  topographic  features  from  multiple  return  LIDAR  data  acquired 
from  airborne  platforms. 

The  project  focused  on  the  first  two  problems  during  the  original  three  year  study,  which  was 
augmented  to  conduct  a  two  year  intensive  study  of  LIDAR  and  hyperspectral  data  for 
coastal  mapping  applications.  The  following  sections  contain  a  review  of  the  approaches 
investigated  and  major  accomplishments  for  the  individual  project  components. 

II.  Research  Results 

A.  Multisensor  topographic  mapping 

Topographic  information  is  provided  by  optical  and  radar  stereo,  single  and  repeat  pass 
interferometry,  and  laser  altimetry.  Stereo-based  optical  methods  are  suited  to  mapping  high 
relief  in  cloud-free  daylight  environments.  Topography  derived  from  stereo  radar  can  be 
acquired  in  a  day-night,  all  weather  environment ,  but  has  poorer  vertical  resolution  than  optical 
stereo.  Interferometry-based  methods  are  superior  in  low  and  medium  relief  areas,  but  repeat 
pass  interferometry  suffers  from  decorrelation  due  to  surface  changes  between  passes,  requires 
precise  information  on  the  position  of  the  sensor,  and  is  affected  by  atmospheric  anomalies. 
Airborne  laser  altimetry  provides  high  resolution  height  information  if  the  position  of  the  aircraft 
is  known  accurately,  but  acquires  data  over  narrow  swaths  that  must  later  be  mosaicked  after 
removing  nonstationary  biases  from  the  individual  flightlines.  New  methods  are  needed  to 
utilize  information  from  DEM’s  acquired  at  multiple  horizontal  and  vertical  resolutions  by  all 
these  methods.  A  capability  to  update  topographic  maps  through  time,  maintaining  the  full 
integrity  and  resolution  of  the  original  data  sets  is  also  required. 

A.l  Dynamic  Adaptive  Grid  Hierarchy.  During  the  first  year  of  the  study,  a  new 
interpolation  method  was  developed  to  integrate  topographic  information  of  different 
resolutions  acquired  from  multiple  sensors  into  a  seamless,  variable  resolution  digital 
elevation  model  DEM,  maintaining  the  highest  resolution  available  for  every  point.  The 
problem  was  formulated  to  maintain  the  integrity  of  the  highest  resolution  DEM  over  an  area 
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and  seamlessly  interpolate  the  DEM’s  at  the  boundaries  using  a  variable  grid  data  structure, 
the  Dynamic  Adaptive  Grid  Hierarchy  (DAGH)  [1,2],  DAGH  was  originally  designed  for 
solving  partial  differential  equations  (PDE)  with  the  capability  of  adaptive  mesh  refinement. 
Whenever  local  errors  from  the  integration  of  the  PDE  on  the  solution  grid  exceed  the 
predetermined  threshold,  it  adaptively  refines  the  area  and  iterates  to  obtain  the  solution 
locally.  The  system  maps  N-dimensional  subspaces  of  different  resolution  into  a  1- 
dimensional  space  using  space-filling  curves  that  achieve  locality  preserving  mappings.  By 
maintaining  records  of  the  relative  coordinates  and  the  instantaneous  resolution  of  subgrids, 
DAGH  then  allows  manipulation  of  data  across  different  levels  while  preserving  resolution 
locally.  The  data  structure  was  modified  for  multiresolution  integration  of  remotely  sensed 
data.  This  required  “inverse”  application  of  the  DAGH  in  the  sense  that  the  locations  of 
refinement  are  predetermined.  Once  a  locality  of  multiple  resolution  data  is  defined,  its 
structure  is  preserved  for  subsequent  analysis.  The  variable  resolution  data  sets  are  then 
represented  with  localized  subgrids  at  different  resolutions,  and  the  Adaptive  Mesh 
Refinement  -  Hierarchical  Data  Format  (AMR-HDF)  is  utilized  to  store  and  retrieve  a  series 
of  multi-dimensional,  multiple-sized,  locally  variable  resolution  data  in  one  file.  The  concept 
is  illustrated  in  Figure  1. 


Figure  1.  Variable  resolution  integration  of  multiple  resolution  information  via  DAGH 

The  methodology  was  implemented  and  evaluated  using  both  simulation  data  and 
multiresolution  data  acquired  over  central  Australia  during  a  NASA  JPL  AIRSAR  campaign 
in  1996.  A  regional  100m  resolution  DEM  based  on  1:100,000  maps,  an  ERS-1  based 
interferometric  DEM,  and  10m  airborne  TOPSAR  data  shown  in  Figure  2  with  the  resulting 
variable  resolution  DEM. 


TOPSAR  (10m)  data  with  fused  result. 
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There  are  several  potential  advantages  of  the  DAGH  approach  for  data  fusion  and/or  variable 
resolution  data  representation:  1)  multiple  data  sets  are  stored  in  a  single  file  with  appropriate 
indexing  using  AMR-HDF;  2)  the  DAGH  allows  potential  reduction  in  storage  when  only  the 
reliable  segments  of  data  are  saved;  3)  by  preserving  the  local  resolution,  seamless  retrieval 
of  data  is  possible;  and  4)  the  ability  to  track  the  evolution  of  local  resolution  enables  the 
extraction  of  uni-grid  (homogeneous  resolution)  DEM’s  as  desired.  The  primary 
disadvantage  of  the  approach  related  to  computational  requirements  for  the  system  and  the 
requirement  for  post-processing  before  ingestion  in  geographic  information  system  software 
packages.  Some  applications  also  require  smoother  transition  across  boundaries  of  data  sets. 
An  alternative  integration  approach  was  developed  to  handle  these  problems. 

A.  2.  Multiscale  Kalman  Smoothins.  During  the  second  year  of  the  study,  the  project  focused 
on  development  of  a  multiscale  interpolation  system  based  on  a  quadtree  framework.  The 
approach  involved  extensions  of  a  traditional  time  domain  state  space  representation 
implemented  as  a  multiscale  model  on  a  pyramid  structure  [5,6,7].  Because  the  models 
satisfy  the  Markov  property  in  scale  and  space,  efficient  recursive  estimation  algorithms 
based  on  extensions  of  Kalman  filtering  and  smoothing  are  employed.  The  covariance  of  the 
errors  is  also  produced  directly  for  a  given  model.  Data  are  represented  on  a  pyramid  of 
multiple  resolutions  as  shown  in  Figure  3. 


Coarse  Resolution 


Figure  3.  Data  pyramid  illustrating  parent/child  relationships  of  MKS  system 

Estimation  proceeds  iteratively.  The  system  is  initialized,  then  upward/downward  sweeps  are 
performed  until  the  estimates  converge.  The  primary  challenges  were  identification  of  the  state 
model,  including  the  transition  matrix  to  maintain  the  Markovian  structure  of  the  model  and 
determination  of  appropriate  measures  of  uncertainty  for  the  model  and  measurement  equations 
that  reflect  the  nonstationarity  of  the  spatial  stochastic  process.  The  methodology  is  fiilly 
described  in  [8,9]. 

The  new  methodology  was  evaluated  using  the  multiresolution  data  acquired  from  central 
Australia,  but  included  3  new  acquisitions  of  5m  TOPSAR  data  over  the  Fink  Gorge  National 
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Park  shown  in  Figure  4.  The  common  area  in  the  TOPSAR  acquisitions  corresponding  to  the 
ERS  subset  shown  on  with  their  corresponding  look  angles.  The  topographic  data  from 
individual  flightlines  is  shown  in  Figure  5  with  the  associated  uncertainties.  Differences  in 
TOPSAR  data  are  due  to  differing  look  angles,  phase  unwrapping  problems,  etc  .  The 
resulting  fused  data  are  shown  in  Figure  6  with  corresponding  estimates  of  local  uncertainties 
and  the  overall  mean  uncertainties.  Figure  7  shows  a  fused  product  over  an  extended  region 
using  ERS  and  TOPSAR  data  at  20,  10,  and  5  meter  resolution  respectively. 


100  km’ 


Figure  4.  Multiscale  regional  topographic  data  over  central  Australia  from  100m  map  based 
DEM,  20m  ERS  interferometry,  and  airborne  TOPSAR.  Look  angles  for  1996  and  2000 
TOPSAR  flightlines  are  indicated. 
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Figure  5.  Subsets  of  2000  TOPSAR  data  from  various  flightlines  with  associated  uncertainty 
acquired  over  the  Finke  River  Ellery  Creek,  and  Palm  Valley  subset  in  2000  (5  m)  resolution. 
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Figure  6.  Fusion  results  and  associated  uncertainties  from  combining  ERS  and  TOPSAR 
flightlines  using  MKS  algorithm 
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Figure  7.  Fusion  results  for  ERS  (20m),  1996  TOPSAR  (10m),  and  2000  TOPSAR  (5m)  data. 

The  methodology  was  also  used  to  fuse  the  30m  DEM  product  derived  from  SRTM  and  5m 
LIDAR  data  acquired  over  Austin,  Texas.  Figure  8  indicates  how  the  method  successfully  filled 
holes  in  the  in  the  LIDAR  DEM  where  vegetation  and  buildings  had  been  removed,  while 
sharpening  the  coarse  resolution  SRTM  DEM. 


and  fused  result. 
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B.  Land  Cover  Mapping  from  High  Dimensional  Data. 

Hyperspectral  instruments  capable  of  acquiring  data  in  hundreds  of  bands  provide  tremendous 
potential  for  discriminating  between  targets  which  are  spectrally  similar.  Data  from  airborne 
systems  were  becoming  more  widely  available,  and  both  military  and  civilian  hyperspectral 
missions  were  scheduled  to  fly  during  the  period  of  the  project.  Standard  supervised  and 
unsupervised  techniques  are  inadequate  for  exploiting  these  data.  During  the  second  year  of  the 
study,  research  was  initiated  on  approaches  for  supervised  analysis  of  hyperspectral  data.  The 
research  focused  on  supervised  classification  of  high  dimensional  input  data  for  problems  where 
outputs  are  also  potentially  high  dimensional  and  the  quantity  of  training  data  is  limited.  A 
multiclassifier  framework,  which  provided  the  capability  to  tune  a  set  of  classifiers  in  binary 
frameworks,  was  adopted  for  the  study.  Two  new  methods,  one  involving  a  pairwise  framework 
and  the  other  utilizing  a  binary  hierarchical  approach,  were  developed.  Within  these 
frameworks,  several  methods  for  band  selection  and  extraction  were  investigated. 

B.  1  Feature  Selection 

High  dimensional  data,  which  are  typically  highly  correlated  within  class  dependent  band 
groups,  are  problematic  for  parametric  classifiers  if  the  quantity  of  training  data  is  limited. 

This  problem  was  addressed  during  the  study  via  various  methods  for  feature  selection  and 
extraction,  which  were  also  implemented  in  combination  within  the  classifier  system 
developed  during  the  project. 

Greedy  forward  selection  was  first  investigated  in  conjunction  with  pairwise  classification 
schemes.  Here,  class  dependent  features  are  selected  based  on  their  incremental  contribution 
to  a  log-likelihood  based  relevance  function  that  provides  a  measure  of  discrimination  of 
each  class  pair  [10].  Because  of  their  capability  to  avoid  local  optima,  heuristic  search 
techniques  were  evaluated  via  an  implementation  of  tabu  search  methodology  [11].  During 
the  last  year  of  the  study,  random  subspace  methods  were  implemented  in  conjunction  with 
bagging  of  training  data  to  address  the  small  sample  problem. 

B.2  Feature  Extraction 

Because  it  maximizes  the  distance  between  populations,  the  Fisher  projection  was  employed  as 
the  primary  feature  extraction  method  for  the  binary  classification  frameworks.  Unstable 
estimates  of  covariance  matrices  were  avoided  by  utilizing  the  extractor  in  conjunction  with 
other  techniques  which  reduced  the  dimension  of  the  input  space  during  a  pre-processing  stage. 
A  new  best  bases  band  combining  technique  which  combines  highly  correlated  adjacent  bands 
for  each  class  was  developed  to  reduce  the  number  of  inputs,  but  retain  interpretable  inputs  [12- 
15].  A  polyline  spectral  approximation  method  was  also  developed,  both  to  reduce  the  input 
dimension  and  provide  a  set  of  stable  features  that  could  be  related  biogeophysical  parameters 
[16]. 
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B.3  Bayesian  Pairwise  Classification 


The  BPC  method  developed  in  this  study  approaches  the  multi-class  problem  with  C(>  2) 


classes  as  a  set  of 


v2y 


simpler  two-class  problems.  For  each  pair  {i,j)  of  classes 


(l  <  /  <  j  <C),  a  separate  classifier  is  trained  to  distinguish  only  between  those  classes. 


Results  of  the  pairwise  classifiers  are  obtained,  and  the  ultimate  class  is  then  selected  either 
via  voting  or  by  the  maximum  Bayesian  posterior  probability  rule  applied  to  an  estimate  of 
posterior  probabilities  obtained  from  the  outputs  of  the  individual  pairwise  classifiers.  While 
the  framework  can  accommodate  a  wide  variety  of  classifiers,  the  approach  was  developed 
with  the  goal  of  successfully  utilizing  weak  classifiers  for  each  pair.  The  system,  which  is 
illustrated  in  Figure  9,  was  implemented  in  conjunction  with  a  Gaussian  model  and  has  been 
used  operationally  for  a  wide  range  of  applications  using  multispectral,  hyperspectral,  and 
radar  remote  sensing  data  [10,  12,  14  ].  One  example  is  included  for  illustrative  purposes. 


final  decision 


Figure  9.  Bayesian  pairwise  classification  framework 


ALI  data  sets  from  2001  (May  31,  June  16,  July  11,  and  August  19)  and  a  single  ETM-i- 
scene  (August  12,  2001)  were  first  classified  as  members  of  the  twenty- three  land  cover 
classes  found  in  the  seasonal  and  occasional  swamps  and  drier  woodlands.  Training  data 
were  selected  manually  using  a  combination  of  GPS  located  vegetation  surveys,  aerial 
photography  from  the  Aquarap  (2000)  project,  and  IKONOS  multispectral  imagery.  Data 
sets  were  randomly  sampled  such  that  50%  were  used  for  training,  with  50%  reserved  for 
testing  the  classifier.  In  order  to  compensate  for  bias  associated  with  co-location  of  training 
and  test  data  within  patches,  additional  sites  were  selected  and  used  for  an  independent 
determination  of  accuracies.  Figure  10  shows  the  resulting  labeled  image  for  the  May  31 
date,  and  Figure  10  contains  classification  results  for  the  scene  using  both  greedy  feature 
selection  and  the  Fisher  extractor.  Detailed  discussion  of  results,  including  comparison  with 
the  standard  maximum  likelihood  classifier,  is  presented  in  [15]. 
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Figure  1 0.  Labeled  ALI  data  over  the  Okavango  Delta,  May  3 1 ,  200 1 
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94.74 

64.34 

85.57 

lloodplaiirl 

67.94 

63.12 

67.94 

61.81 

lloodp[ain2 

61.14 

54.13 

72.54 

62.22 

noodp[aii'i3 

41.71 

34.43 

30.29 

29.78 

water 

99.63 

99.26 

99.63 

99.26 

aL|Lialic  vegclalion 

93.23 

99.44 

95.83 

98.92 

Hrescar2 

95.65 

100 

96.74 

98.89 

Ovfrall 

68.12 

67.87 

knppii 

.66^1 

.663 

Table  1.  Classification  accuracies  for  Okavango  Delta  Data  using  BPC  with  feature  selection 
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B.4.  Bayesian  Hierarchical  Classification 

The  Binary  Hierarchical  Classifier  (BHC),  which  was  originally  developed  external  to  this 
project,  uses  a  Generalized  Associative  Modular  Learning  System  (GAMLS)  algorithm  for 
determining  the  hierarchical  tree  structure  [17,18],  A  hierarchical  tree  structure  is  generated 
by  dividing  the  C  class  problem  into  C-1  two  class  problems  (referred  as  meta-classes),  as 
illustrated  in  Figure  1 1 . 

In  the  original  implementation,  each  meta-class  contains  its  own  unique  feature  space 
extracted  using  the  Fisher(l)  or  Fisher(m)  projection.  The  top  down  BHC  is  initiated  such  that 
the  root  node  contains  all  classes.  Classes  are  allocated  to  two  internal  child  nodes  using  the 
GAMLS  method  [17].  The  process  is  recursively  repeated  as  each  internal  node  is  separated 
into  two  meta-classes  until  each  class  is  fully  separated  as  a  leaf  node.  The  TD-BHC  was 
developed  as  a  decomposition  method  for  classification  problems  with  complex,  high 
dimensional  output  spaces.  It  is  also  advantageous  in  that  it  is  computationally  efficient  and 
less  greedy  than  common  agglomerative  schemes,  including  a  bottom-up  implementation  of 
the  method  (BU-BHC). 


c5-(i . q 


n...-(3i  a.-tci 


Figure  11.  Bayesian  Hierarchical  Classification  Framework 

The  BHC  framework  was  the  foundation  of  the  research  in  hyperspectral  methods  during  the 
last  two  years  of  the  project.  The  small  sample  size  problem  is  critical  for  virtually  all 
hyperspectral  problems  studied.  Sequential  separation  of  the  classes  into  natural  meta-class 
groupings  mitigated  the  impact  of  the  small  training  data  sets  as  the  training  data  for  the 
mixed  classes  was  adequate  for  estimating  parameters  at  higher  levels  of  the  tree.  The  best 
bases  and  polyline  extractors  were  utilized  only  when  the  small  sample  sizes  were 
encountered  at  lower  levels  of  the  tree  and  only  for  those  classes  which  required  reduction  in 
the  input  dimension. 

BHC  has  been  utilized  to  classify  both  airborne  and  space-based  multispectral  and 
hyperspectral  data.  Classification  results  from  corresponding  hyperspectral  scene  acquired 
by  the  Hyperion,  the  hyperspectral  instrument  on  the  EO-1  satellite,  on  May  31,  2001,  are 
shown  in  Figure  12,  and  corresponding  accuracies  for  various  implementations  of  the 
algorithm  are  shown  in  Figure  13.  The  classes  have  been  modified  somewhat  to  reflect  the 
classes  in  the  reduced  area  covered  by  Hyperion  acquisition. 
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Figure  12.  Subset  of  May  31,  2001  Hyperion  data  (bands  51,  149,  31)  and  the  resulting 
classified  image. 


Figure  13.  Sample  mean  and  standard  deviation  of  overall  classification  accuracies  for  BHC 
implementations . 
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C.  High  Resolution  Mapping  from  Light  Detection  and  Ranging  (LIDAR)  Data 

High  resolution  topographic  maps  provide  critical  information  for  a  wide  range  of  scientific, 
public  and  commercial  sector,  and  military  applications.  Capable  of  acquiring  elevations 
with  decimeter-level  accuracy,  LIDAR  is  quickly  becoming  the  technology  of  choice  for 
creating  topographic  maps  over  limited  areas  requiring  high  resolution  information. 

Detailed,  accurate  topography  is  particularly  critical  in  low  gradient  coastal  environments 
which  often  undergo  frequent  change  due  to  environmental  and  anthropogenic  impacts.  In 
the  coastal  zone,  small  changes  in  height  relative  to  sea  level  can  dramatically  affect  local 
geomorphology  and  vegetation  cover,  both  of  which  are  of  interest  to  the  military  for  landing 
in  the  beach  zone  and  navigating  coastal  wetlands.  An  extension  of  the  original  project  was 
funded  by  the  Office  of  Naval  Research  and  conducted  as  a  collaborative  project  with  the  UT 
Bureau  of  Economic  Geology.  This  study  involved  extension  of  an  algorithm  developed  for 
extracting  surface  features  from  multiple  return  LIDAR  and  evaluation  of  its  potential  via  a 
case  study  conducted  over  Matagorda  Island,  a  sandy  barrier  island  on  the  Texas  coast.  The 
science  component  of  the  study  focused  on  obtaining  high  resolution  topography  to  gain  a 
better  understanding  of  the  dynamics  of  habitat  distribution  on  low-lying  barrier  islands. 

C.  1  Feature  extraction  from  airborne  LIDAR. 

The  LIDAR  processing  algorithm,  which  adapts  parameters  to  local  terrain,  involves  a 
sequence  of  operations  shown  in  Figure  14.  Data  points  are  first  filtered  to  remove  extreme 
long  and  short  ranges  that  are  typically  associated  with  anomalies,  and  then  gridded  using  the 
minimum  value  of  the  last  return  in  each  grid  cell.  Low  frequency  variation  in  the  terrain 
represented  by  a  local  average  surface  is  first  subtracted  from  the  pseudo  ground  surface, 
yielding  high-pass  zero-mean  residual  information.  The  ground  values  are  detected  in  this 
“pseudo”  ground  surface  through  a  novel  application  of  the  lower  envelope  follower  (LEE) 
used  to  recover  the  message  signal  from  an  amplitude-modulated  signal.  (In  general,  a  LEE 
tracks  a  modulated  signal  that  decreases  or  slowly  increases.  However,  if  the  modulated 
signal  increases  quickly,  the  LEE  decays  exponentially.)  The  concept  was  adapted  to  the 
problem  of  developing  a  2-D  lower  envelope  (surface)  of  the  LIDAR  data  to  identify 
candidate  ground  level  values.  The  resulting  lower  envelope  surface  is  used  to  threshold  the 
high-pass  residuals,  creating  a  mask  of  the  data  points  in  the  initial  ground  surface  detected  to 
be  on  the  lower  envelope  of  the  signal.  Together,  these  operations  form  the  lower  envelope 
detector.  Unfortunately,  many  ground  points  are  not  detected  in  low  relief  regions  as  the 
LEE  follows  the  lower  envelope  of  the  noise  in  these  areas.  To  recover  these  points,  the 
method  grows  the  lower  envelope  detection  surface  with  a  gradient  flood  fill  algorithm. 

Bridges  and  large  buildings  with  flat  surfaces  may  not  have  been  removed  using  these 
operations,  so  a  gradient-based  feature  detector  is  applied  for  these  man-made  structures.  A 
final  ground  surface  of  the  entire  region  is  created  by  interpolating  values  to  fill  gaps.  This 
surface  is  then  used  to  classify  ground  and  non-ground  points  by  thresholding  the  original 
LIDAR  data  values.  Resulting  non-ground  points  are  separated  into  building  and  non¬ 
building  points  by  detecting  the  points  along  the  buildings’  roofs  and  edges  using  first  and 
last  return  observations.  Details  are  contained  in  [22] 
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Ground  Detection 


Figure  14.  Flow  diagram  for  extracting  bare  earth  topography  from  LIDAR  data 


C.  1.  Matagorda  Island  case  study. 

The  field  area  was  20  km^  of  the  southwest  end  of  Matagorda  Island,  an  undeveloped  barrier 
island  on  the  central  Texas  coast  (Figure  15).  This  area  comprises  an  open-ocean  sandy 
beach,  multiple  dune  lines,  ridge  and  swale  topography,  back  barrier  stabilized  and  active 
dune  fields,  relict  recurved  spits  and  tidal  channels,  and  a  large  relict  washover/flood  tidal 
delta  fan  (Figure  16).  The  area  is  now  part  of  the  Matagorda  National  Wildlife  Refuge, 
although  it  once  served  as  an  air  force  bombing  range.  No  urban  development  has  occurred 
on  the  island  and  the  only  development  present  now  includes  decommissioned  runways,  a 
few  buildings,  and  ditches  and  dikes.  A  DEM  of  a  section  of  the  island  (Figure  17)  was 
derived  and  compared  to  a  map  of  habitats  manually  developed  from  color  infrared  aerial 
photography  and  field  visits  (Figure  18). 


Figure  15.  Location  map,  Matagorda  Island,  Texas 
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Figure  16.  Aerial  photograph  of  the  study  area  looking  northwesterly  toward  San 
Antonio  Bay.  The  surf  of  the  Gulf  of  Mexico  is  at  the  bottom. 


Figure  17.  Shaded  relief  image  of  digital  elevation  model  of  the  study  area.  Banding  in  open- 
water  areas  is  oriented  parallel  to  the  acquisition  flight  lines  and  is  the  result  of  about  0.05  m 
vertical  error  across  the  data  swaths.  The  banding  is  only  apparent  on  the  very  smooth  water 
surface. 


14 


Safi 

AfJtofiio 


Guifot 

M&xico 


□  low  inarsh 
Hlhi^h  marsFi 

□  scnubfthruD 
■  low 

H  his)h  Hal 

Bcoeri’Cicean  beKh 
lemb!>ranly  Haod^d 
_  ?wflle54plowiout5 
■  Iteshwaler  ponds 
IB  upland 

□  ^ann-pamn  onartHy  Nooda'd 
swfliosiblowouts 


Figure  18.  Habitat  classification  map  developed  by  interpretation  of  color  infrared 
photography  and  field  visits. 

Figure  19  is  a  plot  of  the  mean  elevations  and  standard  deviations  of  the  habitats  mapped 
from  the  aerial  photography.  The  units  are  arranged  along  the  x-axis  in  the  expected  order  of 
increasing  elevation,  and  it  is  apparent  that  the  DEM  is  in  accordance  with  this.  Average 
elevations  of  the  intertidal  and  upland  habitats  have  a  total  range  of  less  than  2  m.  Habitat 
elevations  are  separated  in  an  expected  vertical  sequence,  but  standard  deviations  (Figure  17) 
show  overlap  between  environments,  indicating  that  elevation  is  not  the  only  controlling 
factor  on  habitat  type  and  that  vegetation  affects  the  LIDAR  based  elevations.  The  average 
elevation  of  low  marsh  areas  was  only  0.22  m  above  the  water  level  of  ponds  interior  to  the 
relict  flood-tidal  delta  in  the  study  area.  Sedimentation  rates  are  expected  to  be  very  low  in 
these  areas  with  no  open-water  communication  with  the  bay.  Thus,  a  rise  in  relative  sea  level 
of  just  0.22  m  will  expand  the  ponds  and  have  a  profound  effect  on  the  marshes.  Based  on 
sea-level  rise  in  the  bay  since  the  1950’s,  this  amount  of  rise  is  expected  to  occur  during  the 
next  55  years  The  LIDAR  did  not  penetrate  water;  thus  the  sea-grass  elevation  is  actually  the 
surface  of  the  bay  water  during  the  time  of  the  survey.  The  same  is  true  for  the  subtidal 
ponds,  which  are  interior  to  the  washover  fan/flood  tidal  delta  complex  and  have  no  channels 
for  communication  with  the  bay,  and  the  freshwater  ponds  and  flooded  swales  and  blowouts 
interior  to  the  upland  area.  There  is  also  a  lack  of  data  in  places  with  a  smooth  water  surface 
due  to  specular  reflection  away  from  the  aircraft.  At  least  some  of  the  overlap  and  high 
standard  deviations  are  caused  by  vegetation  such  as  the  upland  scrub/shrub  unit.  Greater 
relief  caused  by  ridges  and  swales,  dunes,  and  blowouts,  however,  is  probably  mostly 
responsible  for  the  higher  standard  deviations  in  the  upland  units.  Mean  elevations  gradually 
increase  from  sea  grass  to  high  marsh  and  from  flooded  swales  and  blowouts  to  upland,  but 
an  abrupt  change  of  0.49  m  occurs  between  these  intertidal  and  supratidal  environments. 
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Figure  19.  Average  heights  and  standard  deviations  above  mean  sea  level  for  barrier 
island  habitats. 

Figure  20  contains  a  plot  of  the  MAIOl  ground  transect  (see  Figure  17  for  location).  The 
ground  elevation,  vegetation  height,  and  LIDAR  data  points  that  are  within  1  m  horizontal 
distance  of  the  transect  line  are  shown.  Vertically,  the  LIDAR  data  points  fall  within  the 
vegetation  cover  and  0. 1  to  0.2  m  above  the  ground  elevation  in  low  marsh  areas  and  0.2  to 
0.6  m  above  ground  in  the  upland  area.  The  vertical  scatter  caused  by  the  vegetation  cover  in 
the  low-marsh  area  masks  morphology  with  relief  of  less  than  0.2  m  occurring  across 
horizontal  scales  of  10  m  or  less. 
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Fig.  20.  Ground  surveyed  substrate  and  vegetation  height  relative  to  mean  sea  level.  Lidar 
data  points  within  a  horizontal  distance  of  1  m  from  the  transect  are  also  plotted. 


Delineation  of  the  boundaries  between  low  and  high  marshes  and  between  tidally  influenced 
and  upland  environments  is  a  very  difficult,  but  important  determination.  Analysis  of  the 
LIDAR  based  DEM  shows  that  the  derived  topography  is  accurate  and  detailed  enough  to 
serve  as  an  independent  and  physically  meaningful  layer  in  analysis  involving  other  remote 
sensing  data  such  as  multispectral  or  hyperspectral  imagery.  Vegetation  causes  biases  and 
scatter  in  the  LIDAR  point  data  resulting  in  standard  deviations  that  overlap  between 
environments.  However,  vegetation  is  not  the  only  contributing  factor  to  overlapping 
elevations.  Proximity  to  and  amount  of  communication  with  bay  waters  may  also  control  the 
development  of  low  or  high  marsh,  flats  or  uplands. 
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Reducing  the  scattering  and  bias  effects  of  the  vegetation  on  the  LIDAR  points  would  make 
the  LIDAR  more  useful  in  barrier  island  settings.  As  noted  previously,  the  instrument  used 
for  this  survey  recorded  two  ranges  for  each  outgoing  laser  pulse.  The  resolution  of  these 
measurements,  however,  does  not  allow  distinguishing  ranges  that  are  within  several  meters 
of  each  other.  This  is  useful  in  a  tall  tree  canopy,  but  not  for  low  vegetation  typical  of  the 
study  area.  Future  advances  with  full  waveform  digitization  or  low  photon  return  systems 
may  provide  additional  capability  for  resolving  these  problems. 
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